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Abstract. We discuss some aspects of the search for identities using 
computer algebra and symbolic methods. To keep the discussion as concrete 
as possible, we shall focus on so-called Apery-like formulae for special values 
of the Riemann Zeta function. Many of these results are apparently new, and 
much more work needs to be done before they can be formally proved and 
properly classified. A first step in this direction can be found in [1]. 



1. Introduction 

The Riemann Zeta function is 

oo 1 

(i-i) cw=Ep ®(s)>i. 

k = l 

In view of the "Apery-like" formulae 



00 1 K 00 (_-\\k + l oc 00 1 

(-) c(^3E^, cw-iE^jp ^-r^y 

one is tempted to speculate that there is an analogous formula for C(5), C(6), 
C(7) and so on. The key word here is analogous. For example, extensive 
computation has ruled out the possibility of formulae of the form 

a ™ c ~ i 

C(5)= ^£^(?r c(6)= ^£^(?)' 

where a, b, c, d are moderately sized integers. Such negative results are useful, 
as they tell us it would be a waste of time to search for interesting formulae 
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of a given form. Thus, it would seem there are no corresponding Apery-like 
formulae for higher zeta values. End of story. Consider however, the following 
result of Koecher [2, 3]: 

(L3> c(5)=2 S^(?r" 5 S^(frS?- 

Koecher's formula points up a potential problem with symbolic searching. 
Namely, negative results need to be interpreted carefully, lest they be given 
more weight than they deserve and unnecessarily discourage further investi- 
gation. Also, it becomes clear that symbolic searching is very much limited 
by the need to know fairly precisely the form of what one is searching for in 
advance. 

Koecher's formula (1.3) suggests that one might profit by searching for a 
formula of the form 

^ (_l)fc + l ~ (~l) fc+1 1 (~l) fc + 1 1 

k=l V k I k = l V k ) j = l k=l \ k I j = l 

where ri, r2, r% are rational numbers. The following (conjectured) 1 formula for 
C(7) was found [1] using high precision arithmetic and Maple's integer relations 
algorithms: 

More generally, we have the (conjectured) 2 generating function formula [1] 

fie) y 1 _ 5 y(- 1 -> TT J + zpC 

^ fc 3 (l-z 4 /A: 4 ) 2^ fc 3 (l fc ) l-2 4 /* 4 -_ i 4 --2 4 ' ' 

— 1 k — 1 V k / j — 1 

Note that the constant coefficient in (1.5) gives the formula for £(3) in (1.2). 
The coefficient of z 4 in (1.5) gives (1.4). We arrived at (1.5) by extensive use of 
Maple's lattice algorithms, combined with a good deal of insightful guessing. 
Interestingly, Maple's convert (series, ratpoly) feature played a significant role. 
The reader is referred to [1] for details. 

Comparing our generating function formula (1.5) with Koecher's [3] 

oo 1 oo /_-|\fc+i / j 2 \ k ~ 1 

( L6 ) E £ (1 _ ,2 A 2) = E -j^y {2 + r^T^J gc 1 - ^ 2 /j 2 ) 



1 See note below. 

2 T hese conjectures have subsequently been proved. See Gr anville and Almqvist's 
preprint tittp: / / www. math. uga.edu~andrew/Postscript /BorBrad.pt . 
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raises some interesting issues related to formula redundancy, and which remain 
unresolved. We address certain of these issues in the next section. 

2. Redundancy Relations 

To mitigate the problem of symbol clutter in what follows requires some nota- 
tion. We denote the power sum symmetric functions by 



P r (k) :-- 




Next, we define functions A, \i by 



" ^ (_-\\k+i n 

^-n^) : =EysrII^(*). 



3 = 1 k = l 

n oo 



1 k ) j=l 



fc=l h \k) j = l 



3 = 1 

In the new notation, (1.2) becomes 

r. 9« 

(2.1) C(2) = 3/x(2,P ), C(3) = -A(3,P ), C(4) = ^(4, P ), 

while (1.3) and (1.4) become 

5 5 25 

(2.2) C(5) = 2A(5, P ) - -A(3, P 2 ), C(7) = ^ A (7, ^o) + y A(3, P 4 ), 

respectively. To illustrate the issue of formula redundancy, consider Koecher's 
formula for £(7) [3] which becomes, in our notation, 

(2.3) C(7) = 2A(7,P ) - 2A(5,P 2 ) + |a(3, P|) - ^A(3,P 4 ). 

In view of the second formula in (2.2), the middle two terms of (2.3) must be 
redundant. Indeed, lattice-based reduction shows that 

(2.4) -2A(5,P 2 ) + -A(3, P 2 2 ) = yA(3,P 4 ) + -A(7,P ). 



Although we currently have no real understanding why interrelations between 
A sums such as (2.4) hold, we decided to limit our symbolic search for Zeta 
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function identities in which no such interrelations exist. 3 This was carried out 
by starting with a "full set" of A sums and checking that a relation holds with 
the relevant Zeta value. Now recurse, using the following scheme. From any 
found relation, toss out the Zeta value. If no relation is found amongst the re- 
maining sums, output the relation that held when the Zeta value was included, 
and report it as non-redundant. Otherwise, systematically discard the various 
A sums from the list, until a non-redundant relation remains. Carrying out 
the aforementioned procedure yields the following formulae which evidently 
exhaust the list of non-redundant formulae for each given Zeta value: 

17C(4) - 36^(4, P ) = 5C(4) - 108^(2, P 2 ) = 0, 

7C(6) + 1944^(2, P 4 ) - 1944 M (2, Pf) 
= 215C(6) - 2592^(4, P 2 ) - 3888/x(2,P 4 ) 
= 229C(6) - 2592^(4, P 2 ) - 3888/z(2, P 2 2 ) 
= 1481C(6) - 2592/x(6,P ) - 3888//(2,P 2 2 ) 
= 313C(6) - 648^(6, P ) + 648^(4, P 2 ) 
= 163C(6) - 288^(6, P ) - 432^(2, P 4 ) 
= 0, 

2C(7)-5A(7,P )-25A(3,P 4 ) 
= 4C(7) - 25A(3, P|) + 40A(5, P 2 ) + 225A(3, P 4 ) 
= 22C(7) - 25A(3, P 2 2 ) + 40A(5, P 2 ) - 45A(7, P ) 
= 0, 

72C(9) + 135A(7,P 2 ) - 147A(9, P ) - 60A(5, P|) - 85A(3, P 6 ) + 25A(3, P 2 3 ) 
= 36C(9) - 540A(5, P 4 ) - 96A(9, P ) + 60A(5, P 2 2 ) - 1130A(3, P 6 ) 

+ 675A(3, P 4 P 2 ) - 25A(3, P 2 3 ) 
= 4C(9) + 196A(5, P 4 ) + 32A(7, P 2 ) - 36A(5, P 2 2 ) + 390A(3, P 6 ) 

- 245A(3, P 4 P 2 ) + 15A(3, P 2 3 ) 

= 4C(9) - 20A(5, P 4 ) + 5A(7, P 2 ) - 9A(9, P ) - 45A(3, P 6 ) + 25A(3, P 4 P 2 ) 
= 116C(9) + 68A(5, P 4 ) + 226A(7, P 2 ) - 234A(9, P ) - 108A(5, P 2 2 ) 

- 85A(3, P 4 P 2 ) + 45A(3, P 2 3 ) 
= 0. 



3 Of course, we cannot prove that (1.4) contains no redundancy, since, for example, we 
cannot even prove that £(7) is irrational. 
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No additional formulae other than the formulae given in §1 were found 
for C(2), C(3) and ((5). We discuss additional uniqueness issues in the next 
section. 

3. Uniqueness and £(4n + 3) 

If one extends the list given in the previous section, it becomes apparent that 
£(4n + 3) evidently has a unique representation in terms of A sums of the form 
A(m, P r ) in which r is always a multiple of four. We exploited this observation 
in [1] to arrive at our generating function formula (1.5). Unfortunately, there 
seems to be no sensible selection to make amongst the formulae for Q{An + 1) 
which gives an analogous generating function identity. Our Maple code for pro- 
ducing all possible non-redundant formulae for ^"(13) ran for over two months 
before it was killed. The resulting incomplete file is over three thousand lines 
long and contains hundreds and hundreds of independent formulae. If a gen- 
erating function identity (other than a bisection of Koecher's) for £(4n + 1) 
is found, it is unlikely that it will be discovered by hunting for the appropri- 
ate representatives from the identities for C(9), C(13), etc. and looking for a 
pattern. 

Recall Ramanujan's formulae [4] 

2n + 2 R R 00 jU-4n-3 

A(,^n + d)-W L) (2A . )!(4n + 4 _ 2A . )! 4 Z. e 2 7rfc _ 1 



and 

B2k -E>4n+2-2fc 



2n+l 



1 

2C(4n + 1) = (2vr) 4n+1 — V (-l) fe+1 (2A: - 1)-, , 
sv ' v 1 2n ' y ' (2k)\ (4n + 2-2/c)! 

^e 2 ^-l n^sinh 2 (7rfc)' 

Here, the additional complexity in the 4n + 1 case arises from taking the 
derivative of the appropriate modular transformation formula. Perhaps there 
is an analogous phenomenon operating in the case of these Apery-like identities 
as well. 
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